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一 种 双 水 平 集 模型 分 割 左 心室 膜 的 方法 
李 林 ， 李 军 华 ， 邵 晓 宇 
(南昌 航空 大 学 江西 省 图 像 处 理 与 模式 识别 重点 实验 室 , 南昌 330063) 


摘 要 : 针对 目前 左 心 室内 外 膜 分 割 方法 ， 存 在 分 割 出 的 轮廓 正则 性 差 、 分 割 心室 膜 模糊 边缘 不 完全 和 分 割 效 率 低 
等 问题 。 提 出 新 双 阱 势 函 数 和 各 向 异性 梯度 矢量 流 (AGVF) 改进 水 平 集 模型 能 量 函 数 ， 用 融合 0 水 平 集 入 水 平 
集 的 双 水 平 集 ， 同 时 分 割 左 心室 内 外 膜 的 方法 。 首 先 ， 用 改进 Hough 变化 圆 检测 算法 ， 定 位 心室 内 外 膜 初 始 位 置 ; 
然后 利用 双 水 平 集 模 型 同时 分 割 内 外 膜 。 观 察 和 对 比分 析 实 验 结果 ， 该 方法 能 够 分 割 出 平滑 的 左 心室 内 外 膜 轮廓 ; 
能 够 分 割 出 的 符合 临床 定义 的 左 心 室内 外 膜 轮 廊 ， 且 分 割 MRI 图 像 左 心室 内 外 膜 轮 廓 重合 率 平 均 提 高 到 0.9569。 
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Double level set model for segmentation of left ventricular membrane 


Li Lin, Li Junhua, Shao Xiaoyu 
(Key Laboratory of Image Processing & Pattern Recognition in Jiangxi Province, Nanchang Hangkong University, 
Nanchang 330063, China) 


Abstract: In view of the existing endocardium and epicardium contour of left ventricle segmentation methods, there are 
some problems such as the poor regularity of segmented contour, segmenting the fuzzy edge of ventricular membrane 
incompletely and low segmentation efficiency. This paper proposes a new double well potential function and anisotropy 
gradient vector flow (AGVF) energy function improved level set model, with the fusion of zero level set and k level set of 
the double level set, the division of left ventricular inner and outer membrane method at the same time . First of all, with the 
modified Hough change detection algorithm, positioning ventricular endocardium and epicardium initial position; Then, the 
double level set model is used to segment the endocardium and epicardium simultaneously. The experimental results show 
that the proposed method of this paper can segment smooth endocardium and epicardium contour of left ventricle. The 
analysis of experimental datum demonstrate that the method of this paper can segment endocardium and epicardium contour 
of left ventricle confirms to the clinical definition, And increases averagely the overlap rate of the endocardium and 
epicardium contour on MRI images to 0.9569. 
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Lo 0 ”引言 灰 度 分 布 差 别 大 等 问题 ， 导 致 这 些 分 割 算 法 分 割 左 心室 内 乡 
二 二 膜 模糊 边缘 时 ， 会 出 现 提取 边缘 不 符合 临床 要 求 且 正 则 性 较 


近年 ， 心 血管 疾病 的 高 死亡 率 ， 致 使 心脏 左 心室 内 外 膜 差 等 问题 。 
分 割 成 为 研究 的 热点 。 常 用 方法 包含 闵 值 法 由、 聚 类 法 叫 、 针对 分 割 心室 内 外 膜 弱 边 缘 正 则 性 差 ， 小 梁 肌 和 乳头 肌 
区 域 生 长 法 中 、 主 动 轮廓 模型 法 所、 水 平 集 法 外 和 动态 规划 方 的 分 割 不 彻底 问题 。 提 出 用 融合 改进 Hough 圆 检测 定位 的 双 
法 四。Lee 等 人 D 中 利用 区 域 生 长 和 迭代 阔 值 法 分 割 做 心室 内 水 平 集 模 型 方法 分 割 左 心室 内 外 膜 。 用 改进 Hough 圆 检测 进 
膜 ， 用 提取 的 内 膜 边界 引导 主动 轮廓 模型 提取 外 膜 ;， 李 晓 守 行 定位 ， 定 位 的 轮廓 作为 内 外 膜 分 割 初始 位 置 ， 用 融合 各 向 
等 人 [融合 Mean-Shift 聚 类 算法 的 视觉 显著 性 方法 对 心 内 膜 异性 梯度 矢量 流 (anisotropic gradient vector flow，AGVF) 
初 定 位 ， 利 用 GVF Snake 模型 分 割 内 外 膜 ， 朱 错 等 人 外用 结 的 双 水 平 集 模型 分 割 出 内 外 膜 轮 万。 

合 K 均值 的 自 适应 核 函 数 带 宽 的 Mean-Shift 算法 分 割 图 像 ， RE i 

区 域 生长 法 提取 目标 区 域 ，Dharanibai 等 人 0 利用 动态 规划 。 “1 ”心室 腊 自 动 分 割 万 法 

分 割 心室 外 膜 用 融合 Otsu 算法 的 主动 轮廓 模型 分 割 心室 内 1.1 心室 膜 自 动 分 割 流程 

膜 , ， 王 兴 家 等 人 00 提 出 用 典 入 先 验 知识 的 耦合 水 平 集 算法 本 文 提出 自动 分 割 左 心室 内 外 膜 方法 主要 步骤 : 包含 改 
模型 分 割 左 心室 内 外 膜 。 由 于 人 体 心 脏 的 复杂 性 ， 在 临床 上 进 Hough 变换 定位 和 双 水 平 集 模型 分 割 , 分 割 方法 流程 图 如 
定义 肌 小 粱 和 乳头 肌 属 于 左 心室 内 腔 ， 且 认为 只 要 腔 室 凸 出 图 1 所 示 。 为 简化 心室 内 外 膜 分 割 流程 ， 解 决 人 为 干预 手动 
部 分 的 周 长 超 过 50% 被 心肌 包围 ， 就 认定 为 该 部 分 心室 为 心 设置 模型 分 割 初 始 位 置 问题 ,采用 改进 Hough 变换 圆 检测 方 
肌 的 一 部 分 ， 每 个 心跳 周期 影像 数据 量 巨大 ,容易 出 现 伪 影 、 法 ， 定 位 心室 内 外 模型 分 割 初始 位 置 ， 针 对 心室 膜 模 糊 边 缘 
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和 冲 陷 区 域 ， 提 出 融合 新 双 阱 势 (new double-well potential， 传统 Hough 方法 是 对 检测 出 边缘 像素 点 进行 遍历 运算 ， 
NDWP) 和 AGVEF 的 双 水 平 集 模型 方法 分 割 心室 内 外 膜 。 最 后 输出 累加 像素 最 多 的 圆 。 图 2 中 ，(Oj 是 未 输入 圆心 大 至 
立 置 ， 对 边缘 检测 图 遍历 运算 后 输出 的 定位 结果 ，(d) 是 输入 


输入 图 像 ”一 一 pp Hough 定 位 和 | 双 水 平 集 模型 


心 位 置 ， 得 到 的 左 心室 内 腔 定位 结果 。 这 种 遍历 运算 耗 时 
较 久 ， 且 检测 结果 定位 效果 较 差 。 
传统 Hough 变换 圆 检测 自 适 应 定位 , 定位 结果 不 满足 要 
求 ;， 输入 圆心 位 置 后 对 内 腔 进行 定位 ， 依 然 无 法 避免 人 为 对 
收敛 性 检 分 割 心室 膜 初 始 位 置 的 干预 。 为 解决 人 为 干预 圆心 种 子 点 设 
置 ， 提 高 定位 精度 和 效率 。 改 进 Hough 变换 圆 检测 方法 ,在 


六 


图 1 心室 内 外 膜 分 制 流程 图 对 图 像 进行 找 圆 定位 前 ， 对 输入 图 像 进行 OTSU 自 适 应 阀 值 
Fig.1 Theflow chart of the endocardium and epicardium 分 割 处 理 ， 像 素 较 大 点 连 成 区 域 设 为 前 景区 域 ， 累 加 器 先 累 
segmentation 加 前 景区 域 的 像素 点 ; 为 避免 检测 出 多 个 无 用 圆 ， 利 用 图 像 
1.2 定位 模型 初始 轮廓 边缘 信息 进行 约束 ， 在 有 限时 间 间 隔 内 只 允许 沿 着 梯度 方向 
传统 水 平 集 模 型 需要 人 为 设置 分 割 初始 位 置 ， 造 成 分 割 的 像素 点 进行 投票 031。 改 进 后 半径 沿 着 边缘 点 (x,y) 法 线 方向 
流程 复杂 且 低 效 , 容易 受 人 为 主观 性 影响 , 为 克服 这 一 问题 ， 变化 ， 而 不 是 所 有 方向 的 遍历 变化 。 圆 心计 算 公式 为 
提出 用 图 像 边缘 信息 约束 的 改进 Hough 变换 圆 检 测 方法 自动 J 0) 
获取 初始 轮廓。 b=y-rcos0 
传统 Hough 变换 0 的 算法 原理 是 利用 点 与 线 的 的 对 偶 es 
性 ， 将 图 像 空间 的 线条 变 为 空间 的 聚焦 点 ， 从 而 检测 给 定 图 。 站/O) ，Q 0 是 对 应 方向 的 榜 度 。 计 算计 
像 是 否 存在 给 定性 质 的 曲线 。 圆 的 方程 为 程 中 ， 
(x-a) +(x-a) = (1) in- 
攻 中 。 (ea 为 圆心 ，* 为 半径 。 首 先 用 筑 子 对 图 像 进行 | ee G) 
中 : (a,5) 为 圆心 ，7 为 半径 。 首 先 用 算 子 对 图 像 进行 边缘 cosO = G 
检测 得 到 二 值 图 像 ，Hough 变换 会 将 图 像 空间 对 应 到 (a.b.7) es 
参数 空间 ， 取 和 图 像 平面 一 样 的 参数 平面 ， 以 二 值 图像 中 每 改进 Hough 变换 圆 检测 的 具体 步骤 : 
一 个 非 零点 为 圆心 ， 用 设 定 的 半径 在 参数 平面 上 画 圆 ， 符 合 a) 根据 左 心室 图 像 特征 , 利用 先 验 知识 确定 半径 的 变化 
(4,b,7) 的 点 进行 累加 ， 即 对 圆 覆 盖 的 坐标 点 进行 投票 ， 随 着 区 间 ; 
半径 步 长 一 直 循 环 运算 到 设 定 半径 的 最 大 值 ， 最 后 找 出 参数 b) 对 输入 图 像 进 行 OTSU 自 适 应 阔 值 分 割 处 理 ， 提 取 
平面 上 的 峰值 , 得 到 (a,b,7) 值 及 其 在 图 像 对 应 位 置 显示 形状 ， 出 像素 较 大 点 连 成 的 前 景区 域 ; 
而 后 通过 调整 六 值得 到 目标 圆 形 。 传统 Hough 变换 中 有 三 个 c) 利用 图 像 边缘 信息 , 约束 前 景区 域 像素 点 沿 着 边缘 梯 
自由 参数 (a.b,r) ， 三 个 参数 在 三 维 空间 中 进行 投票 累加 ( 设 度 方向 进行 投票 ; 
自 适 应 圆 检测 半径 范围 为 R )， 则 总 共 进 行 了 RN; 投票 计算 ， d) 输出 所 有 像素 投票 后 ,， 累加 数组 中 值 最 大 像素 对 应 的 
复杂 度 函 数 为 : f(N)=0(N) 。 其 中 ，N 为 图 像 边缘 检测 后 ， 圆 。 
边缘 像素 点 数 。 对 改进 后 Hough 变化 圆 检测 方法 进行 复杂 度 分 析 。 由 于 
5 SA 图 像 边 缘 信息 的 约束 ， 沿 着 梯度 方向 的 像素 点 进行 投票 ， 变 
成 了 三 个 自由 参数 在 二 维 空间 中 进行 投票 ， 则 共 进 行 了 RN? 


投票 计算 ， 算 法 复杂 度 函 数 为 : fi(N)=O(N?) 。 其 中 ， N 为 
图 像 经 过 OTSU 自 适 应 阔 值 分 割 处 理 后 ， 像 素 为 1 的 像素 点 
总 数 。 

对 比 图 2 和 3 可 以 看 出 ,使 用 改进 后 方法 定位 效果 较 好 ， 
定位 结果 被 整个 内 膜 包 围 ， 提 高 了 对 内 膜 的 定位 精度 。 对 比 


j(X) 和 f(N) ,改进 后 检测 圆 算法 ,降低 了 算法 空间 复杂 度 ， 
(a) 原 图 (b) 边 缘 检 测 结 果 提高 了 运算 效率 。 原 Hough 变换 圆 检测 耗 时 3.209 秒 ， 改 进 


(a)Original image (b)Edge detection results 后 只 需要 0.3920 秒 就 能 准确 定位 到 心室 内 膜 。 


(c) 自 适应 定位 (d) 输 入 圆心 后 定位 (a)Otsu 处 理 (b) 定 位 结果 
(c) Adaptive localization (d)Location after input center (a)Otsu processing (b)Locating result 
图 2 Hough 变换 圆 检测 妈 3 ”优化 算法 定位 过 程 


Fig.2 Hough transform round detection Fig.3 Localization process of optimization algorithm 
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1.3” 双 水 平 集 模型 分 割 方法 
1.3.1 水 平 集 算 法 原理 


水 平 集 (1 的 算法 思想 就 是 用 高 维度 表示 低 维 度 ， 在 高 维 
看 上 闭合 曲线 轮廓 。 在 传统 的 


面 用 隐 函 数 等 值 点 来 表示 平 
平 集 算法 模型 中 ， 出 现 反 复 初始 化 、 曲 线 结果 正则 性 差 等 
题 ，Li 等 人 [5 提出 了 一 种 距离 正则 项 水 平 集 (distance 
regularized level set evolution，DRLSE) 模型 算法 ,解决 了 传 


统 水 平 集 模型 反复 初始 化 问题 。 根 据 曲线 的 运动 规律 ， 定 义 
能 量 泛 函 为 
E($)=1R,($)+ ec () (4) 
其 中 : R,($) 是 定义 的 距离 正则 项 ，4 是 大 于 零 的 常数 ， 外 着 
能 量 项 so.(y) 由 图 像 先 验 信息 或 图 像 底层 视觉 特征 决定 。 
水 平 集 的 距离 正则 项 被 定义 为 
B,(#)=| pVoDar (5) 


其 中 : 7 为 势 函 数 ， 
potential，DWP) 定义 为 


D:[0,oo) >R。 


传统 势 函 数 (double-well 


EAA 
nllveD=1( (6) 
JVI-1), Ivok1. 


偏 微分 方程 为 
[ae Ivk1 

$ 

ot 


or (Ca) IvE1 (2) 


扩散 率 为 
sinc(2|VYg) |Y 吕 < 和 1 
nv vg (8) 
|VYg| 


1.3.2 新 双 阱 势 函 数 
水 平 集 函数 的 正则 性 在 水 平 集 曲 线 演变 中 是 内 在 保持 
的 ， DRLSE 水 平 集 分 割 模糊 边界 得 到 的 边界 正则 性 太 差 ,为 
得 到 更 加 平滑 的 分 割 轮廓 曲线 ,在 水 平 集 模型 中 提出 NDWP 
代 蔡 传统 势 函数 约束 模型 曲线 演化 。 
定义 新 双 阱 势 函数 NDWP 为 


34vwldvgl-DP， voi 
| 9) 
30voP -D-iogdygl) [VoP1. 


pa(|V#))= 


如 此 定义 ， 知 双 阱 势 函数 有 两 个 最 小 值 分 别 是 |Y4F0 和 
|YgFI1。 


NDWP 扩散 率 为 
(21V¢$|-D(UVGI-D |Y 刘 <1 

户 (|VY 人 |) = gp "ene Ivdlz1 (10) 
偏 微分 方程 为 

8 [VI(21V¢$| -DUVGI-D) IvVAs!1 

2 11) 

(i Mr logioe |Y 外 >1 ( 
且 oo = (12) 


证 明了 势 函 数 的 扩散 速率 是 有 界 的 。 得 到 一 个 NDWP 
的 距离 正则 化 项 ， 水 平 集 的 演化 是 一 个 梯度 流 ， 它 可 以 使 能 
量 函 数 最 小 化 ， 且 对 水 平 集 函 数 的 梯度 模 进 行 约束 ， 避 免 传 
统 水 平 集 曲线 演化 中 的 重新 初始 化 。 由 图 4 可 得 

a) 当 |Yg| >1， 线 正 向 扩散 ， 越 扩散 1Y4| 越 


»>n>0,， 


等 : 一 种 双 水 平 集 模型 分 


ChinaXiv 合 作 期 刊 


害 左 心室 膜 的 方法 


小 ， 曲 线 运 动 到 |Ygy>1， 
DWP 模型 方法 ; 
b) 当 1/2<Vg| <1， 
越 大 ，NDWP 模型 方法 分 割 效率 高 于 DWP 分 割 方法 ; 
co) 当 |Vb| <1/2 ， 曲 线 正 向 扩散 ， 越 扩散 1Y4| 越 小 ， 
演化 的 的 结果 是 1YY%F0，1Yg| <1/2 区 域内 距 模型 方法 科 
线 较 近 ， 模 型 扩散 速率 对 曲线 演化 影响 较 小 ; 
总 之 ， 曲 线 扩散 向 着 使 |Ygh>0 或 1Y4h>I 的 方向 进 
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NDWP 需要 的 先 


i<n<0, 
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失 代 次 数 明显 少 于 


线 反 向 扩散 , 越 扩 散 1vy| 


线 
等 值 


一 一 号 
4 一 一 op 


3 a 
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(a) 双 阱 势 (b) 扩散 速率 
(a) Double well potential (b) Diffusion rate 


图 4 双 阱 势 函数 及 对 应 扩散 速率 


Fig.4 Double well potential function and corresponding diffusion rate 


1.3.2 构造 AGVF 模型 
于 DRLSE 模型 在 分 割 凹陷 区 域 时 ， 


容易 出 现 从 分 割 


现象 ， 
分 割 


很 难 正 确 分 割 出 目标 轮廓 。 因 梯度 矢量 流 09 能 将 模型 
线 演 化 引 向 边缘 深 目 处， 为 使 模型 有 较 好 分 割 效果 ， 


在 水 平 集 模型 中 提出 用 AGVF 改变 图 像 的 外 力 场 分 布 。 


构造 扩散 矩阵 ， 定义 结构 四 通道 张 量 Eh 


Td | u? uu, 
| T_ | “x “ry 
1， J Vu Vu 
bE [UU Us 


道 进行 非 线 性 处 理 : 
Oy . i 
cet Dv vw jv | 


大 ,1=1 


~ 


对 各 个 通 


ovuP)= 05] 


分 解 非 线性 结构 张 量 ， 由 于 / 是 正定 对 称 阵 ， 
(7-4E)x=0， 其 中 为 单位 矩阵 ，x* 是 非 零 向 量 ， 得 至 
和 切线 方向 对 应 特征 值 4 。 


1 
4 = tt (tJa) +4 启 ] 
1 
= +4 局 ) 
再 利用 入 入 ， 


am =0 
(1 -hE)1=0 


求 相应 特征 向 


如 


1 Mu+Vo -hb 
i Ja 本 


人 即 


LK 
[2 J ~ J» 中 


根据 特征 值 对 其 进 


加 


N= 
ge a9) 
Ci+(I-C)exp(-C:/( 和 一 入) 和 关 思 


其 中 Cs(0.D,C >0 常数 。 定 义 扩散 矩阵 为 


D= NH + ho po 


则 推出 各 向 异性 扩散 方程 为 


(13) 


(14) 


(15) 


利用 
1 法 线 


(16) 


(17) 


(18) 


(20) 
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史 =4div(Dvu) (21) 
同 理 推出 : 
Sdn(DV) (22) 


根据 图 像 信 息 构建 AGVF 模型 ， 设 1(%y) 是 左 心室 内 外 
膜 边 缘 梯 度 ， 左 心室 图 像 梯度 矢量 流 为 
V(%y)=(u(y)v(%y)) ， 对 其 分 量 进行 处 理 ， 得 到 极 小 化 能 量 
泛 函数 0 


e=|| uate 


其 中 :oo 为 wwY 对 *》 的 一 阶 偏 导 数 : 4 是 由 图 像 噪声 
控制 的 参数 ， 与 噪声 成 正比 关系 。 通 过 求解 变 分 法 推出 的 欧 
拉 方 程 组 ， 可 以 得 到 GVF 场 ; 
Hu-(u-f Rt 
| -A(R+A 
其 中 : 六 是 拉 普 拉 斯 算 子 ,控制 边缘 梯度 量 向 周围 空间 扩散 ; 
(u- 太 ) 和 (vw- 太 ) 保 真 项 , 主要 为 保护 图 像 每 点 像素 梯度 值 保 


w+ HV Pv -VP drdy (23) 


(24) 


持原 图 一 致 ， (万 + 万 ) 为 保 真 系数 ， 在 分 割 区 域 边界 处 取得 
最 大 值 。 给 定 一 点 (%y)， 由 欧 拉 方程 组 推出 能 量 方 程 为 


e=||s (VDGa+t t+)+ 
Qi-g(v (Cf) (uA) )ardy 


其 中 : g(|Vf|)=exp(-(|VA/K))，K 梯度 敏感 参数 。 
GVF 的 数值 实现 UL9 知 ， 当 a 时 ， 保 证 GVF 能 


量 泛 函 是 收敛 的 。 其 中 At 是 时 间 步 长 ， AcAy 为 空间 步 长 。 
可 以 看 出 ， 当 图 像 比较 粗糙 〈 即 Ax 或 Ay 较 大 ) 和 控制 参数 4 
较 小 时 ， 收 敛 速度 较 快 。 
为 了 减 小 AcAy 对 收敛 速度 的 影响 ， 在 GVF 模型 中 引入 
图 像 的 各 向 异性 ， 在 GVF 能 量 泛 函 中 加 入 各 向 异性 函数 
div(DVu) 和 div(DVv) ， 为 更 好 的 收敛 到 图 像 边 缘 凹 陷 处 ， 添 
加 约束 分 割 区域 和 边缘 的 函数 gs 和 h(x) 


(25) 
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AGVF 处 理 图 像 ， 扩 大 了 模型 的 扑 捉 


(b) 原始 外 力 场 分 布 


(b) Distribution of the initial external force field 


(a)U 型 图 
(a)U shape image 


同居 


9 


(d)AGVF 外 力 场 分 布 
(d)Distribution of the AGVF 


external force field 


(c)GVEF 外 力 场 分 布 
(c)Distribution of the GVF 
external force field 
图 5 外 力 场 分 布 
Fig.5 Distribution of external force field 


1.3.3 双 水 平 集 能 量 函 数 


由 于 单一 水 平 集 曲 线 演 化 模型 分 割 心室 内 外 膜 边 缘 ， 需 
要 进行 两 次 分 割 ; 且 ! 于 心室 内 外 膜 边 缘 像素 点 对 比 度 个 同 ， 
同一 模型 分 割 内 外 膜 ， 内 外 膜 会 出 现 边缘 泄漏 。 为 提高 模型 
分 割 效率 和 模糊 边缘 分 割 精度 ， 提 出 用 改进 后 水 平 集 模型 的 
0 水 平 集 和 K 水 平 集 融 合 的 双 水 平 集 模型 ， 同 时 分 割 心室 内 
外 膜 。 
双 水 平 集 方法 的 能 量 函 数 为 

&.($)=LUR(9)+AL(9)+aA(9y) (28) 
其 中 : R(9) 模型 距离 正则 项 , L(y) 是 模型 边缘 约束 项 , A(9) 是 


g(|VvU|)=exp(-(|AU;/K))， h(|V0;|))=1-g(|VU;|)，K 梯度 敏感 
参数 ， 天 值 与 图 像 梯度 相同 处 的 收敛 速度 相关 ， 天 值 越 大 ， 
梯度 值 相同 处 的 收敛 速度 就 越 快 ，U; 分割 区 域内 第 i 类 信息 
和 矩阵。 得 AGVF 模型 能 量 泛 函 ; 


E=||(gdvywD)+sg(YFD)(div(CDY+div(CDY)) 


(26) 
+(h|UiD)+ pun(| VADGV -VU;) dxdy 


Un Un 人 ye 
7=| 2 让 | 根据 上 述 扩散 箱 隆 构造 过 程 推出 p 。 令 


U=U;， 当 AGVF 能 量 函 数 取 极 小 值 时 ， 推 导出 AGVF 模型 
的 计算 方法 为 


P=(s (lvU D+ sllv DaivDvu)- 
(avo D+ p(w ud,) 
=(s(lvUD+ slv I) (div(Dv)- 
(alvo D+ en Dv-0,) 
通过 以 上 公式 推导 ,理论 上 本 文 提出 的 AGVF 模型 处 理 
图 像 时 扑 捉 范围 更 大 。 图 5 用 U 型 图 来 验证 提出 模型 的 优越 
性 。 图 5(b) 图 像 原始 外 立场 凹陷 处 呈 水 平 状态 ， 通 过 比较 ， 
图 5(c) GVF 外 力 场 改 进 外 力 场 范围 ， 但 没有 图 5(d) AGVP 
外 立场 的 凹陷 区 域 收 剑 性 更 强 ， 范 围 更 大 。 对 比 发 现 引入 


(27) 


模型 演化 速度 控制 项 。R(9) 保持 双 水 平 集 模型 演化 过 程 中 | 
线 的 正则 性 ; L(y) 融合 AGVF 约束 0 水 平 集 和 水 平 集 朝 着 
内 外 膜 边缘 及 边缘 凹陷 处 演化 ;4(p) 约束 模型 曲线 朝 目标 边 
缘 速 率 。 定 义 R(o) ，Z(p) ，4(O) 如 下 : 


R(O)=|pdVgl)dx (29) 


L(p)=4|g5. (Gp) voldxt hl 825.((0) -Rv¢ldx 30) 


A($)=a |gH,(-p)dxt+o|gH,(k-p) dx GD 
且 5 和 分 别 为 狄 拉克 函数 和 赫 维 赛 德 函数 : 
. COS 4 
-| [' (2 ls G32) 
0, Xl>e. 
1 rx,l1 7 
2 E+rLsin[ 到 】 ， xlze 
和 H(zx)=11, X>e (33) 
0, 4 < 一 人 
边缘 指定 符 为 
2 
SGrlvo iF 
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G; 为 标准 偏差 o 的 高 


域 。 


有 f 核 函数 ，1 的 定义 在 整个 图 像 区  ” 理 过 程 中 , 本 文采 
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模型 分 割 方法 ， 


实验 参数 | 


双 水 平 集 模型 将 水 平 集 演化 导出 为 梯度 流 ， 用 优化 距离 KH=02，, Ah4=5, m=-30，, %h=6, 
正则 化 项 和 外 部 能 量 来 最 小 化 能 量 泛 函 ， 从 而 驱动 零 水 平 集 
向 目标 边缘 运动 。 距 离 正 则 项 约束 拓扑 水 平 集 演化 具有 唯 件 平台 为 Windows7， matlab2014a。 
的 向 后 和 向 前 扩散 效应 ， 在 水 平 集 附近 有 符号 距离 分 布 , 它 ”2.1 分 割 结果 
能 够 维持 水 平 集 函 数 的 期 望 形状 ， 约 束 模型 对 模糊 边缘 的 分 本 文 分 别 用 CT 图 像 和 MRI 
市 。 距 离 正则 化 消除 了 传统 水 平 集中 需要 的 重新 初始 化 ， 简 。 例 一 为 CT 图 像 分 各 


司 时 保证 应 有 的 计算 精度 。 


引入 AGVF， 促 使 模型 完全 分 割 凹陷 边界 ， 减 少 迭 代 次 数 的 ”(a) 是 左 心室 原 图 
室内 外 膜 结果 ，(c) 是 利用 双 水 平生 


改进 Hough 变换 圆 检 测定 位 和 双 水 平 集 


多 次 实验 所 得 ， 设 置 ==15 ， 
=-40 。 实 验 配置 处 理 


器 为 IntelXeon(R)@1.60GHz 1.60 GHz，RAM 8 GB 开发 软 


图 像 展示 分 割 效果 ， 其 中 案 


| 结果， 案例 二 是 MRI 图 像 分 割 结果 。 用 
此 了 运算 。 选 用 有 限 差分 法 来 量化 图 像 梯度 ， 外 部 能 量 函数 MRI 图 像 阐 述 本 文 方法 分 割 


心脏 内 外 膜 主要 过 程 。 图 6 


双 水 平 集 模型 分 割 方法 主要 部 分 的 伪 代 码 如 下 ; 


效果 ，(d) 是 分 割 


获得 定位 内 外 膜 的 初始 轮廓 像素 点 数据 集 ; 
计算 出 图 像 像 素 点 到 内 外 膜 初始 轮廓 的 欧 氏 距离 ; 


While《〈 双 水 平 集 能 量 泛 函 取得 最 小 值 ) % 人 迭代 nm 次 


{ 


距离 正则 项 函数 约束 初始 轮 局 


曲线 沿 法 线 方向 扩散 ， 


位 半径 为 Rr ， 则 双 水 平 集 获 得 的 初始 像素 点 数 共 有 


图 像 边缘 梯度 控制 曲线 演化 速度 ; 
} a 
对 双 水 平 集 模 型 算法 进行 复杂 度 分 析 。 设 内 外 膜 初 始 定 (a) 原 


1=2x(R+r) ， 设 图 像 像 素 点 数 为 m， 则 计算 图 像 像素 点 到 初 


始 轮廓 欧 氏 距 离 的 时 间 应 为 =/(m-/) 。 进 入 迭代 循环 内 ,1 
线 演化 过 程 中 选用 有 限 差 分 方式 进行 表示 ， 


(9 凶 一 各 )/At=L(t)， 其 中 At 是 时 间 步 长 ，ZL( 如 ) 为 曲线 第 n 


次 演化 的 像素 点 集合 。 曲 线 每 次 演化 过 程 中 涉及 到 像素 点 数 
为 (2 9 线 迭 代 过 程 中 的 时 间 复杂 度 TD = 钦 / 十 狼 / 十 … 侈 》 左 


心室 内 外 膜 分 割 结果 推出 ， 


<n* 加 。 因 为 如 <m， 所 以 双 


水 平 集 模型 分 割 的 时 间 复 杂 度 为 : T(m)=0(m) 。 (c) 分 割 


水 平 集 方法 就 是 求解 


个 曲线 随时 间 变 化 的 偏 微分 方 


程 。 而 双 水 平 集 模型 是 在 三 维 空间 上 ， 找 到 满足 0 水 平 集 和 


k 水 平 集 曲线 点 的 集合 ， 


所 以 双 水 平 集 的 空间 复杂 度 为 ，5(d)=0(4?) ，4 代表 维 数 。 


2 ”实验 结果 与 分 析 


0 将 三 维 问题 利用 二 维 知识 解决 。 


效果 的 三 维 


图 


(a)Original image 


效果 


(c)Segmentation effect 


Sb 于 
能 量 图 。 


(b) 内 外 膜 定位 图 


改进 Hough 变换 圆 检测 定位 心 
站 模型 方法 同时 分 制 内 儿 


膜 


(b)Location map of endocardium and 


epicardium 


onhsnm oo 


所 印 


150 100 50 50 100 150 


(d) 分 割 结果 能 量 图 


(d)Energy diagram of 


the segmentation results 


图 6 本 文 方法 分 割 过 程 


Fig.6 The segmentation process of this method 


为 了 更 好 地 证 明 本 文 算法 的 精准 、 高 效 性 ， 做 了 多 组 对 


比 实验 ， 即 与 经 


RSF 模型 


本 文 实验 数据 集 采 


日 了 和 鲁 昂 大 学 计算 机 工程 与 科学 。 轮廓 模型 (Snake)、 


ja | 


MRI 心脏 影像 库 P5 和 ”中 原 大 学 CT 心脏 影像 库 ， 和 鲁 昂 大  ” 结果 进行 对 比 。 


学 数据 库 包 含 48 人 心脏 序 


图 7(a) 和 


多 


每 个 心脏 振动 周期 有 387 引 


20、CV 模型 、 基 于 区 域 的 主动 
传统 DRLSE 模型 方法 分 割 左 心室 内 外 膜 
8(a) 手 动 分 割 结果 作为 标准 。 

列 , 每 份 包含 279 张 图 像 , 图 像 大 图 7(b) 和 图 8(b) 是 Snake 模型 方法 分 割 结果 ;图 7(c) 和 图 8(c) 
小 为 256x216 ;中原 大 学 数据 库 包 含 一 心脏 的 10 个 振动 周期 ， 是 DRLSE 模型 方法 分 割 结果 ; 名 


8(d) 是 RSF 模型 分 割 方法 


(a) 手 动 分 割 


(a)Manual segmentation 


(b)Snake 模型 半 


长 图 像 ， 图 像 大 小 256x256 。 程序 处 结果 ; 图 8(e) 是 CV 模型 方法 分 割 


结果 。 
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(b)Semi-automatic segmentation of Snake model 


(c)DRLSE 模型 
(c)DRLSE model 


(d) 文 献 [21] 
(d)Document [21] 


(e) 本 文 方法 
(e) Method of this paper. 
图 7 案例 1 分 割 效果 图 


Fig.7 Segmentation effect of the case one 


(a) 手 动 分 割 (b)Snake 模型 分 害 


(a)Manual Segmentation (b)Snake model segmentation 


Nt 


(c)DRLSE 模型 分 割 (d)RSF 模型 分 割 
(c)DRLSE model segmentation (d)RSF model segmentation 


(e)CV 模型 分 割 中 本 文 方法 
(e)CV model segmentation (e)Method of this paper 


图 8 案例 2 分割 结果 


Fig.8 Segmentation results of case two 
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2.2 实验 结果 分 析 表 2 案例 2 分 割 结果 对 比 
为 定量 分 析 模 型 分 割 左 心 室内 外 膜 效 果 ， 引 入 重 琶 率 Table 2 ”Compared segmentation results of case two 
( Dice ) 和 误 分 割 率 〈 ME ) 2 量化 分 析 实 验 结果 ， 客 观 对 方法 模型 Dice 
比 各 模型 方法 分 割 效果 ; Snake 模型 0.9244 
Die 多 | | (35) DRLSE 模型 0.9438 
| Fanoa | + | Fro RSF 模型 0.9298 
ME 1 [Bra NB |+ | Fmt Funw| (36) CV 模型 0.9125 
| Byuna |+ | Frma | 本 文 方法 0.9569 
其 中 :Re 和 Bu 是 专家 手动 分 割 左 心室 外 膜 的 目标 区 域 其 次 ， 为 验证 本 文 算法 具有 良好 的 分 割 效 果 ， 分 割 MRI 
和 背景 的 像素 集 ， Fs。 和 Bw 是 用 方法 分 割 左 心室 外 膜 的 目 图 像 时 , 加 入 了 基于 边缘 分 割 的 RSF 模型 方法 和 基于 区 域 的 
示 区 域 和 背景 的 像素 集 ，|.| 集 合 基 数 表示 某 一 集合 所 属 的 。” CV 模型 方法 作对 比 。RSF 模型 方法 分 割 边缘 时 ， 模 糊 边 缘 
类 。 约束 条 件 较 弱 ， 出 现 边缘 过 分 割 现象 ，CV 模型 方法 分 割 模 
直接 从 视觉 角度 分 析 图 7 案例 1。 基 于 区 域 的 Snake 模 糊 区 域 时 ， 出 现 多 个 拓扑 结构 ， 边 缘 发 生 泄漏 。 由 图 8 直接 
型 方法 分 割 对 比 明显 边界 时 ， 分 割 效果 恨 好， 分割 模糊 边界 “从 视觉 角度 就 能 看 出 RSFE、CV、Snake、DRLSE 模型 方法 分 
误差 较 大 ; 模糊 边缘 需要 人 工 二 次 分 割 校正 ， 图 7(b) 较 粗 蓝 。 制 左 心室 外 膜 时 ,图 像 像 素 间 对 比 度 明 最 区 域 效 果 良 好 , 但 在 
线 为 手动 校正 二 次 分 割 位 置 ， 虽 然 Li 引入 势能 函数 解决 传 心室 膜 没 有 明确 界定 处 ， 分 割 所 得 边缘 正则 性 差 ， 甚 至 在 弱 
统 水 平 集 模型 演化 过 程 中 反复 初始 化 问题 ， 但 模型 对 初始 位 ”边缘 处 出 现 边缘 泄露 现象 。 对 比 图 8(f) 本 文 算法 分 割 的 左 心 
置 比较 敏感 且 仍 需要 手动 设置 ， 分 割 弱 边缘 凹陷 处 时 ， 易 出 “室外 膜 结 果 ， 本 文 算法 由 于 引入 AGVF， 在 边缘 凹陷 区 域 也 
现 欠 分 割 现象 ; 文献 [21] 方 法 在 模糊 边缘 处 出 现 过 分 割 现象 ; 行 充分 分 割 , 新 双 阱 势 函 数 NDWP 的 提出 , 使 距离 正则 
本 文 方法 分 割 边缘 为 类 圆 形 符合 左 心室 外 膜 临 床 定 义 ， 没 有 项 能 够 更 好 的 约束 边缘 正则 性 ， 得 到 符合 临床 定义 的 平滑 类 
ee 边缘 泄露 现象 。 避 形 轮廓 ; 比较 表 2 多 种 模型 方法 分 割 左 心室 , 得 到 的 Pice 数 
© 表 1 不 同方 法 分 割 左 心室 据 ， 发 现 本 文 提出 方法 分 割 左 心室 内 外 膜 Dice 高 至 0.9565， 
© Table ] Segmentation of left ventricle by different methods 得 到 了 高 精度 的 左 心室 内 外 膜 轮 万。 
© 算法 程序 运行 时 间 Dicee ME 通过 上 述 分 析 ， 本 文 算法 虽然 空间 复杂 度 任 和 水 平 集 模 
© 40.0521 0.9478 0.0261 型 是 同一 数量 级 ， 但 在 分 割 时 间 上 ， 于 是 对 内 外 膜 进 行 同 
i Snake 模型 半自动 分 害 400364 09513 00379 时 分 割 ， 耗 时 仅 约 为 单 水 平 集 分 割 的 一 半 ， 自 动 初始 定位 ， 
©O 39.2635 0.9486 ”0.0279 降低 了 模型 对 初始 位 置 的 敏感 性 ， 自 动 分 割 左 心室 内 外 膜 ， 
O) 39.2973 0.9457 0.0158 提高 了 心室 内 外 膜 弱 边 缘 处 和 贴 陷 区 域 的 分 割 精度 ， 人 避免 J 
SS 20.7844 0.9001 0.0261 人 工分 割 的 主观 性 影响 .图 7 和 8 证 实 本 文 算法 无 论 分 割 CT 
QQ 30.4342 0.8933 ”0.0249 心室 图 像 还 是 MRI 心室 图 像 , 都 能 有 效 的 完成 左 心室 内 外 膜 
DRLSE 模型 ee 
a 26.2163 0.9054 0.0308 的 自动 分 割 。 
:一 27.8563 0.9038 ”0.0302 3 ”结束 语 
SC 33.1733 0.9175 0.0348 
5 32.6621 0.9235 0.0441 本 文 提 出 基于 双 水 平 集 模型 的 左 心室 内 外 膜 的 分 割 方 
人 文献 2] 31.2753 0.9246 0.0371 法 , 利用 NDWP 提高 模型 分 割 左 心室 膜 的 演化 速度 , 平滑 模 
31.2315 0.9241 0.0361 糊 边 缘分 割 轮廓 ; 用 加 入 AGVF 的 外 部 能 量 函 数 驱 动 曲 线 引 
© 15.2223 0.9337 0.0198 向 凹陷 区 域 边 缘 ， 自 动 分 割 左 心室 内 外 膜 。 根 据 心脏 左 心室 
本 17.6135 0.9334 0.0201 医学 影像 的 特点 对 图 像 进行 定位 , 利用 改进 Hough 变换 检测 
本 文 方法 14.8946 0.9355 0.0168 到 的 圆 形 ， 作 为 双 水 平 集 模型 的 初始 轮廓 ， 通 过 双 水 平 集 模 
15.6873 0.9362 ”0.0173 型 分 割 得 到 左 心室 内 外 膜 轮 廓 。 与 Snake 模型 算法 、RSF 模 
结合 表 1 案例 一 数据 ， 本 文 提 出 方法 用 二 维 平面 解决 三 。 型 算法 和 CV 模型 算法 的 分 割 结果 相 比 ， 本 文 算法 无 论 分 割 
维 空间 问题 , 虽然 模型 空间 复杂 度 与 DRLSE 模型 和 文献 [21] ”CT 图 像 还 是 MRI 图 像 ， 保 证 分 割 的 内 外 膜 边缘 的 正则 性 ， 
模型 仍 属 于 同一 数量 级 ， 但 在 分 割 时 间 上 明显 少 于 DRLSE ” 均 提 高 了 内 外 膜 的 分 割 精度 和 分 制 效率 ， 人 避免 了 初始 轮廓 的 
模型 和 文献 [21] 的 分 割 时 间 ， 大 幅度 提高 了 分 割 效率 ; 在 分 。 人 为 设置 ， 实 现 了 左 心室 外 膜 的 自动 化 分 割 。 然 而 ， 观 察 分 
制 精度 方面 ，Snake 模型 半自动 手动 二 次 校正 分 割 ， 分 制 结 ”市 过 程 发 现 ， 内 膜 边缘 像素 间 对 比 明显， 应 充分 利用 局 部 信 
果 Dice 平 均 才 只 有 0.9484， ME 平均 为 0.0269; DRLSE 模型 ” 息 优化 内 膜 分 割 。 接 下 来 工作 针对 模型 不 足 进行 优化 ， 提 高 
方法 分 割 结 果 Dice 平 均 为 0.9006, ME 平均 为 0.0280; 文献 [21] ” 本文 方法 分 割 精度 ， 构 建 左 心室 三 维 模型 。 
分 割 结果 结果 Dice 平 均 为 0.9224，ME 平 均 为 0.0380; 本 文 参考 文献 : 
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